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We present preliminary results concerning amplitude saturation in (3 Cephei 
Q-i models. Using nonlinear approach we have investigated amplitude limitation 
q mechanism in (3 Cephei stars. In our approach radial modes have been treated as 
!— i ' representative for all acoustic oscillations. We have studied pulsation properties 
c/2 . of several models (7-20 M Q , Z = 0.02, 0.015) using radiative Lagrangean 
' hydrocodes (essentially those of Stellingwerf 1974, 1975). Nonlinear limit cycles 
, (monoperiodic full-amplitude oscillations) have been calculated through Stel- 
lingwerf s (1974) relaxation technique, which also provides information about 
limit cycle stability. 

In our main survey (Z = 0.02) only the fundamental and the first overtone 
modes are linearly excited. Nonlinear growth rates have been used to determine 
the modal selection (see e.g. Stellingwerf 1975). We found that fundamental 
mode pulsation is dominant. First overtone pulsation is restricted to inter- 
mediate masses and to the vicinity of the blue edge. The first overtone and 
the fundamental mode pulsation domains are separated by either-or or narrow 
double-mode domains. 

Predicted, single mode saturation amplitudes have been compared to am- 
plitudes observed for monoperiodic (3 Cephei variables (Fig. 1, left). Predicted 
amplitudes are significantly higher. The amplitudes may be lowered by decreas- 
ing the metal abundance of the models, Z. We have found that for Z = 0.015, 
decrease of model amplitudes is not sufficient. At the same time instability strip 
shrinks and leaves a lot of stars beyond the blue edge. By lowering Z we are 
not able to match simultaneously the observed amplitudes and the instability 
strip. 

Predicted amplitudes may be easily lowered to the observed level if one 
assumes collective saturation of the pulsation instability, by n similar acous- 
tic modes. In this hypothetical multimode solution, amplitudes of individual 
modes are a factor of ~ y/n lower than in single mode solution. Using linear 
code of Dziembowski (1977) we have determined the number of linearly unsta- 
ble acoustic modes for models of different masses, located in the center of the 
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Figure 1: Predicted single mode saturation amplitudes (left) and amplitudes assuming 
collective saturation by several acoustic modes (right). For comparison, amplitudes 
of monoperiodic (open circles) and multiperiodic (filled circles) j3 Cephei stars are 
ploted. 



main sequence band. This number doesn't vary much along evolutionary track 
and thus, was assumed to be representative for all models of a given mass. The 
number of unstable modes is much higher than the number of detected modes 
in the multiperiodic (3 Cephei variables. Nonlinear simulations (Nowakowski 
2005) also show that not all unstable modes take part in the saturation pro- 
cess. Thus, we have arbitrarily assumed that only one third of linearly unstable 
modes take part in the saturation. Amplitudes rescaled under assumption of 
collective saturation are presented in Fig. 1 (right). Using only part of the lin- 
early unstable acoustic modes, we have lowered theoretical amplitudes to the 
observed level. Thus, we argue that collective instability saturation is sufficient 
to explain observed amplitudes of the (3 Cephei pulsators. A possible difficulty 
of this model is that predicted pulsation-induced broadening of spectral lines 
might be higher than observed. We discuss this problem in Smolec & Moskalik 
(2006). 

In several of our radiative models we have found numerically robust, double- 
mode behaviour, with radial fundamental and first overtone modes simultane- 
ously excited. This form of pulsation is encountered only in intermediate mass 
models (10-11 M Q ). Depending on the specific model, the origin of double 
mode pulsation can be traced to one of two different mechanisms: either to 
the nonresonant coupling of the two excited modes, or to the 2w\ ~ ujq + u>2 
parametric resonance. 

Full results of this analysis (including discussion of non-uniform filling of 
the theoretical instability strip by j3 Cephei variables and detailed study of the 
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double-mode models) are presented in Smolec & Moskalik (2006). 

Acknowledgments. The authors are grateful to the EC for the estab- 
lishment of the European Helio- and Asteroseismology Network HELAS, which 
made the participation of the authors at this workshop possible. This work has 
been supported by the Polish MNil Grant No. 1 P03D Oil 30. 

References 

Dziembowski, W.A. 1977, AcA 27,95 

Nowakowski, R. 2005, PhD Thesis, Copernicus Astronomical Center, Warsaw 
Smolec, R., Moskalik, P. 2006, MNRAS, submitted 
Stellingwerf, R.F. 1974, ApJ 192, 139 
Stellingwerf, R.F. 1975, ApJ 195, 441 



